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Abstract 

We consider the numerical evaluation of one dimens i onal p rojections of general multivariate 
stable densities introduced by lAbdul-Hamid and NolanI l)l998j) . In their approach higher order 
derivatives of one dimensional densities are used, which seem to be cumbersome in practice. Fur- 
thermore there are some difficulties for even dimensions. In order to overcome these difficulties 
we obtain the explicit finite-interval integral representation of one dimensional projections for all 
dimensions. For this purpose we utilize the imaginary part of complex integration, whose real part 
corresponds to the derivative of the one dimensional inversion formula. We also give summaries on 
relations between various parameterizations of stable multivariate density and its one dimensional 
projection. 



1 Introduction 

Stable distributions are known as the limiting distributions of the general central limit theorem and 
the time invariant distributions for Levy processes. Because of these probabilistically very impor- 
tant properties, many studi es have been carr i ed ou t on their theoretical as p ects. Comprehensive 
books have been published ((Arauio and Ginel ()l98nl ) or IChristoph and (jl993l ^ ) . Furthermore 
applications of stable distributions to heavy-tailed data have been growing over past few decades. 
These applications have appeared in many fields like finance, Internet traffic or physics. Especially 
in time series, the observations from the causal stable ARM A or frac tional stable ARIMA model ca n 
be cons i dered as a multivariate sta ble random vector. We refer to Uchaikin and Zolotarev ( 19991 ) . 
Rachev (j20n,'^ or Adler et al. (j2nfl,'j ) for these applications. As various multivariate heavy tailed data 
have become available in these fields, the treatments of multivariate stable distributions are needed. 
However, the multivariate stable distributions require handling of their spectral measures which is 
computationally difficult. 

We briefly review some important recent progresses concerning the multivariat e stable distribu- 
tions. Concerning the calculation of the integral with respect to the spectral measure. iBvczkowski et al 
have approximated the spectral measure by a discrete spectral measure which is uniformly close 
to the original spectral measure. The y also have given several inst ructive graphs of the stable densi- 
ties. Concerning the approximations, Davvdov and NagaevI ( 2002l l have outlined possible directions 
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althou gh t heir results are still theoretical. F or the estimation of the parameter s, see Pivato and Secol 
( 20031 ^ and Davvdov and PaulauskasI ( 19991 ^ and for the hypothesis testing see Mittnik et al.l ( 1999^ . 
Although there are many other researches, accurate calculations of the multivariate stable densities 
are clearly of basic importance. 

Fo r the calculations of the multivariate stable densities, the method of Abdul-Hamid and NolanI 
199S ) based on one dimensional projection is very promising. In the approach of Abdul-Hamid and NolanI 



19981 ^. the basic ingredient is a function ga,d{v,(3) defined by 



1 ro< 
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1 ro< 
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[(3 tan 



u 



d~i. 



du, 



(^vu + ^(3uloguj u'^~^e~'^du, a = 1, 



(1) 



where a is the characteristic exponent of the stable distribution, d is the dimension and /? corresponds 
to the skewness parameter. 

Evaluating ga,d{u,f3) using this very definition is not very practical because the integrand is in- 
finitely oscillating and changing its sign. Therefore an alternative evaluation of ga,div, (3) is desirable. 

In this paper, we obtain explicit finite-interval integral representations of ga,div,P), which do 
not require the derivatives of one dimensional densities. The integrand in our representation is well- 
behaved without infinite oscillations. Furthermore, our integral representation covers all dimensions 
and all parameter values, except for an exceptional case of a < l,/3 = 1 and even dimension, where 
we need an additional term which is easily tractable. 



Th e calculation of ga,d{v,P) corresponds to the evaluation of the equation (2.2.18) of IZolotarev 
|l98()l ). which is stated without an explicit proof. For odd dimension, the real part of the equa- 
tion is equivalent to ga,d{v,(3) and for even dir nension, the imag inary part is equivalent to ga,d{v,P). 



Complete derivations of the equation (2.2.18) of lZolotarevI ()198(tI ) is also useful for the numerical evalu- 
ations of the derivatives of one dimensional symmetric densities investigated in Matsui and Takemural 
(20Qi). 

This paper is organized as follows. F or th e rest o f this section we prepare preliminary results 
based on lAbdTil-Hamid and Nola^ (Il998h and iNolanI (Il998h . In Section HI we derive finite-interval 
integral representations of one dimensional projections. In SectionlHlwe summarize relations between 
various parameterizations of the multivariate stable density and its one dimensional projection as 
propositions. Some discussions and directions for further studies are given in Section^ Some proofs 
are postponed to Section [3 



1.1 Definitions and prefiminary results 

Let X = {Xi, X2, ■ ■ ■ , Xd) denote a d-dimensional a-stable random vector with the characteristic 
function 

^>„,d(t) = £;exp(i(t,x)), 

where t = {ti,t2,... ,td) and ( , ) denotes the usual inner product. Among several definitions of 
stable distributions, we give the spectral representation of ^a,d(i)- This requires integration over 
the unit sphere S'^ = |s : l [s|| = 1} in M*^ with res pect t o a sp ectral measure P (Theorem 2.3.1 of 
Samorodnitskv and Tannul |l994l V Theorem 14.10 of lSatnl (|l 9991 ^1. The spectral definition is given as 
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2) 



'exp (- /g, |(t,s)|° (1 -isign((t,s))tan^) T{ds) + i{t, fi)) , a / 1, 
^exp (- /g, |(t, s)| (1 + if sign((t, s)) log |(t, s)|) T{ds) + i(t, /i)) , a = 1. 

The pair (F, /i) is unique. The definition ((2j correspo nds to Zolotarev's (A) representation and we 
use the notation Saj{T,fi) for this definition following NolanI (|l998l ^. 

For another definition, we give a multivariate version of Zolotarev's (M) parameterization, which 
is also defined in iNolanI (|l998l V We use notation S^^{r,^o). 



exp(-/g,|(t,s)|"(l + isign((t,s))tan^(|(t,s)|i-"-l))F(ds)+i(t,^o)) 
exp (- /g, I (t, s) I (1 + if sign((t, s)) log | (t, s) |) F(ds) + i{t, //q)) , 



a 7^ 1, 
a = 1. 



(3) 



In the equation © , we avoid the discontinuity at a = 1 and the divergence of the mode as a — > 1 which 
are observed in the definition in the non-symmetric case. The difference of (M) parameterization 
from (A) representation is only 



Mo = M 



vra f „ / , N 
tan — / s F(as), 

2 J§d 



Now consider inverting t he cha r acteri stic f unction ^n^t) to evalua t e the multivariate stable 
density. The useful idea of NolanI ( 19981 ^ and Abdul-Hamid and NolanI ( 19981 ) is to perform this 
inversion with respect to the polar coordinates. The projection of an a-stable random vector onto 
one dimens ional stable random variab l e is e xplained as follows. Here we follow the arguments of 
Ex. 2. 3. 4 of ISamorodnitskv and Tannul (|l994l l. For any t G R*^ consider the characteristic function 
E exp{iu{t,x)) of a linear combination Y = '^■tiXi, i.e., 



E exp{iu{t, X 
^ fexp (- 
I exp (- 



/^)) 

" [/g, |(t,s)rr((is) -isign(tx)tan^/g,sign((t,s))|(t,s)rF(ds)]) , a ^ I, 
[/g, I (t, s) \r{ds) + if sign(tx) /g, (t, s) (log + log I (t, s) I) r{ds)] ) , a = l. 



Then, from the definition below, ii^ exp(iii(t, x)) is considered as a characteristic function of one 
dimensional stable distribution with the following parameters. 



a(t) 

/3(t) 
^(t) 



a{t)) 
0, 



\{t,s)rT{ds)y\ 

sign(t,s) |(t,s) 



T{ds) 



a / 1, 

f /g,(t,s)logKt,s)| r(ds), a = l. 



(4) 
(5) 
(6) 



Regarding u as the length and t as the angle, we transform the d-dimens ional rectangular integral into 
the polar coordinate integral. In the following theorem of iNolanl ()l998l V multivariate stable densities 
of (A) representat ion are projec ted onto one dimensional densities in (A) representation. Since in the 
L'iginal theorem dNolanl (|l998l )) there are some trivial typographical errors, we correct them in the 



or 



formula Q and the following theorem. 
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Theorem 1.1 (jNolanI (|l998h ) Let < a < 2 and X = {Xi,X2, . . . ,Xd) be an a-stable random 
vector Sa^di^: ^) with d>2. Then the density fa,di^) ofJ^ is given as 

,/3(s)V(t(s))-'^ ds, a^l, 




a(s) 

where ga,d{v,P) is given in ^\). 

Abdul-Hamid and NolanI ( 19981 ^ expressed ga,d{v,P) as an integral over a finite interval, utilizing 



the integral expressions of one din iensional densitie s. They explained the method when d is odd 
as follows. From Theorem 2.2.3 of IZolotarevI (|l986l ^. the one dimensional density ga,iiv,P) can be 
written as ^ 

ga,i{v,P)= hg{9;a,P,v)d9, 

J a 

where —tt/2 < a < h < it /2 and kg is a somewhat complicated but explicit function. Then, for odd 
d, differentiating g^^i d — \ times with respect to v gives the representation 

gaAv^P) = c{d) j Ij^^^y dO, 

where c{d) is some constant. However, the computation of d'^~^hg{9;a, f3,v)/dv'^~^ seems to be 
cumbersome. Furthermore, they did not give the corresponding expression for even d. In the follow- 
ing, we obtain explicit finite-interval integral representations of ga,d{v,(3), which do not require the 
derivatives of one dimensional densities. 



2 Finite- interval integral representations of projected one dimen- 
sional functions 

In this section we obtain the finite-interval integral representations of ga,d{v, (3) in It is based on 
the function h^[x]a, j3) defined below in which corresponds to the n-th order derivative of one 
dimensional inversion formula including the imaginary part. 



We first prepare some notations which are the same as in IZolotarevI (|l986l l. 

K{a) = a — 1 + sign(l — a), a / 1, 
e = (3K{a)/a, a ^ 1, 

^ sin f + e) Y/('~») COS f {(g - l)y + aO} 
Ua{^;9) = [ , a /I, 

\ cos ^(p J COS ^(p 

0) = 2 ^ exp [-^[^ + ^) tan j , 

r /sin»(,+../2)y/a-) 

I exp(-x//3 + ((^ + 7r/(2/?))tan93), a = 1, 
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(a/x)V(i-"), a^l, 
exp(— X — 1), a = 1, 



Vni^) = r" [^if) {r (|(/?) sin ^(n + 1)((^ + 1) + r (^(/.) cos ^(n + 1)((^ + 1)} 
= r" (^(/.) {r (|(^) sin |(n + l)(c^ + 1) - r' (^(/.) cos ^(n + 1)((^ + 1)| 



We now define a function a, /?) as 



TT 

where 



/i"(x;a, /3) = — I (iz)" exp {izx + iIj{z; a, —(3)) dz, (7) 







V'(2;a,/3) 



-z" exp(— z^6'a), a 7^ 1, 
— log z, a = 1, 



is the logarithm of the characteristic function of stable distributions which are analytically extended 
from the semi-axis Re z > into the complex plane. Then, ga,d{v,(3) in is expressed in terms of 
a, 13) as follows. Note that j3 in ga,d{v, (3) is not equal to (3 in h^^x; a, f3), because /? in h'^{x; a, j3) 
is based on one dimensional (B) parameterization. Therefore, we need the following definitions, 

Pb = — — — — ■ arctan ( /3 tan — I, 6b = l3BK{a)/a, a^l. 



nKia) V 2 



Proposition 2.1 Define x = [cos{'k /2a9B))^^°' v and y = tt/2v + /31og(7r/2). Then a representation 
of ga,div, (3) in ^ in terms of h{- ;a,/3) is given as follows. ■ 



ga,d{'"^P) 



Re ^h''~\y*;a,\P\), 



a 



1, 



where (3*^ = Pb signx and y* = y sign/3. 

Proof of this proposition is given in Section |SJ Although the imaginary number i is present in 
the above formulas, it is for notational convenience and we do not have to write Re/i(-) or lmh{-) 
depending on the dimension d. 

By Proposition 12. 11 it suffices to express h^(x; a, (3) as a finite-interval i ntegral. We now give our 
main theorem. We mention that the following results a,re partially state d in Zolotarev ( 19861 ) for the 
real case without proofs and the equation (2.2.34) of IZolotarev (jl98fj ) is not complete. Therefore 
we give the complete results including the calculation of imaginary parts. The full proof is given in 
Section 



Theorem 2.1 Let x > in the case a 7^ 1 and (3 > in the case a 
integral representations of h"'{x;a, f3) are as follows. 
(a) (a / 1, /3 / 1) or (a > 1, /3 = 1) 

"1 



1. Then the finite-interval 



Re 



exp 



-x'^/^'^-^^Ua{ip;9)) VMd^. 
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Im h^{x;a,p) = ^ j'^exp (-x"/("-^)C/,(v9; 0)) Wn{ip)dip. 

(b) (a = l,/?y^l) 

Reh^{x;l,P) = ^y"'^exp(-e-"/%i(<^;/3)) K(<^)dc^. 
lmh^{x;l,p) = ^y'^exp(-e-"/%i(v9;/3)) Ty„((^)(iy;. 

(c) (a<l,/3 = l) 

Re h^{x\a, 1) is the same as in (a). 

Im /i"(x;a,l) = -/" exp f-a;"/(°-^)C/^(99; 0) W„(¥;)d(^ - - / exp(xr - r°)r"dr. 

(d) (a = !,/? = !) 

Re 1, 1) is the same as in (b). 

Im 1, 1) = exp (|-e~^/^f7i(v3; W„((/j)(ivJ - - ^ exp (xr + r log r) r"(ir. 

Interestingly, th e case /3 = 1 is different from other values of the parameters. As suggested by 



ZolotarevI ()198(tI ) , the representations of Theorem 12.11 are essentially different from the n-th order 



deriva tive of the integral representations of densities given in the equation (2.2.18) of IZolotarev 

3 Relations between representations of multivariate stable densities 
and one dimensional projections 

In this section, we present propositions concerning the projections of a stable random vector onto 
one dimensional stable random variable. The idea of projection is useful for the purpose of the 
computation of the multivariate densities. There exist two representations (A) and (M) for the 
multivariate stable distributions, whereas three representations (A), (B) and (M) are known for 
general one dimensional stable distributions. The projection of the multivariate repr esentat i on fA ) 
onto one dimens i onal r epresentation (A) is given in Theorem ll.il which is obtained bv lNolanI ( 19981 '). 



Although iNola 3 (!l998l ^ mentioned the projection of (M) onto (M) as a theorem without proof, his 



is different from ours. Therefore we describe the projection of the multivariate representation (M) 
onto the projected one dimen sional represe ntation (M) with a precise proof. Note that the idea of 



projection entirely belongs to iNolanI ()l998l l. Further we present results for all projections without 
proofs by adding the other 4 projections. Note that for a = 1 multivariate representations (A) and 
(M) and one dimensional representations (A) and (M) are the same. Therefore we omit statements 
on the projections (A ^ M), (M ^ A) and (M ^ B) for a = 1. 
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3.1 One-dimensional projections in various representations 

Here we present propositions on various projections. Note that the parameters cr{t), /3(t) and /u(t) 
used in this section have aheady been defined by Q), © and © in Section [TTTl respectively. 

Proposition 3.1 (The projection (M M)) Let < a < 2 and X = {Xi, X2, . . . , X^) be an 

a-stable random vector S^^^{T,i') with d>2. Then, the density f^^^i'^) is given as follows. 

/5W = //M ( ^""''47''^'^ /3(s)) (a(s))-'^ ds, 



where 



9aA^,P) = (/3tan^(n-^.-))) u'-'e-^'^du 

fiZit)=t^n^(^m^it)- ljt,s) rids 



and 



When a = 1, (A) parametrization in Theorem II . II and (M) parametrization in Proposition 13. II are 
the same. It is easy to confirm /j'^(x) = /i.d(x) and gf^^{v,f3) = gi.d{v,(3) only by notational change. 
Note that at a = 1, /ijf (t) becomes 

/^M(t) = -f/3(tMt)loga(t)- / (t,s) log|(t,s)|r(ds)' 
= /z(s) + (2/7r)/3(s)cT(s)loga(s). 



Proposition 3.2 (The projection (A — > B)) Under the same conditions and notations as in The- 
orem Yn\. the density fa,di^) is given as follows. 



/a,d(x) 



' Is^9^4 {M^)r'' ds, a / 1, 

Is^9f,, ( <--;;-)-y^(-) ,/3(s)) {asis))-' ds, a = 1, 



where 



vu — u 



(2^) 



Kd /o°° cos {vu + f3u log n) n"' ^e ''/^"du, 



a = 1, 



and 



aB(t) = { (cos(f X(a)fe(t))) 
|-(t), 



— arctan ( /3(t) tan 
7rK{a) V 2 



17^, "7^1, 
a = 1, 



/ii?(t) =C7B(t)/3(t)logaB(t)-- /" (t,s) log|(t,s)|r((is) 
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Proposition 3.3 (The projection (A M)) Under the same conditions and relations as in The- 
orem [7771 the density fa,d{^) is given as follows. 

where 

TTOi 

//A/(t) = o-(t)/?(t) tan — . 



Proposition 3.4 (The projection (M A)) Under the same conditions and relations as in Propo- 
sition rOl the density /(^^(x) is given as follows. 

fU^) = [ 9.4 f ^"""''lT^^^'\ /?(s)) (a(s))-'^ ds, a / 1, 
//^(t) = -tan^ / (t,s) T{ds). 



Proposition 3.5 (The projection (M B)) Under the same conditions and relations as in Propo- 
sition WH the density /*^^(x) is given as follows. 

/.'f.(x) = / g^a f ^"""''^",^^'^'\ fe(s)) Ms))-' ds, a + 1. 



3.2 Finite-integral representations for other representations 



In IZolotarevI (|l986l ) only the analytic extension of the characteristic function defined by (B) repre- 
sentation is considered and Theorem 12.11 corresponds to (B) representation of stable distributions. 
Therefore, we have introduced the parameters /3b and Ob in the finite integral-representation for (A) 
of Proportion 12. 11 Alt hough i t is po ssible to derive (A) or (M) representation versions of Theorem l2.1l 
utilizing the results in iNolanI \\mi^ . we have to consider the analytic extension of the characteristic 
function defined by (A) or (M) representation which needs very complicated arguments. Therefore, 
in this paper, we confine our results to the case of (B) representation and for (A) and (M) represen- 
tation, we utilize Proposition 12.11 and I.S.HL Note that, accordingly, for g^di'^^f^) Proposition l',^.7\ 
we do not need to use the extra parameters like Pb or 6b other than that defined in Theorem 12. II 

Proposition 3.6 Let g^^ ^{v , fi) he as in Provosition \3.1\ Definex = (cos (7r/2a6'B))^^" (f -|- tan (7r/2Q;^B)) 
and y as in Provosition \2.1\ Then the representation of g^^^{v, f3) using function h{- ■,a,/3) is as fol- 
lows. 



9^A^,P) 
where (3*b = 13 sign 2: . 



{ Re (-(:/;f/_r^" /^'^-n|x|;a,/3|,), a / 1, 

a = 1, 



2<^(7ri) 
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Proposition 3.7 Let g^^ ^(v , /3) be as in Proposition XS.^ Then the representation of g^^{v,l3) using 
function h{- ; a, (3) is as follows. 

where (5* = (3 sign?; and v* = v sign/?. 



Re 



Re 



]^h<'-^{v*;l,\/3\), a = l, 



4 Some discussions and future works 

In this paper we focused on one dimensional projections g of the general multivariate stable density. 
This is only one step in calculating the density itself. We need to substitute parameters which are 
functions on the unit sphere s € S*^ into g and then we need to integrate over the unit sphere to 
evaluate the density itself. We presented improvements in evaluation of g. Improvements of other 
steps are also important for the evaluation of the multivariate density. By further careful examinations 
of the integrand of the function h^{x; a, (3), we migh t find some us eful regularities like one dimensional 
finite integral representations stated in Section 3 of 

For our future work we consider showing various properties like tail dependencies or relations 
between the spectral meas ure and tai l s. The ore tically the asvm ptotic estimates of multivariate stable 
densities are obtained by IWatanahel l^2Q0(i ] or iHirahal (|2q3). Relations between these theoretical 
results and our representations or numerical results are of our next concern. 

Furthermore, we can also consider the expansions or the asvmpt otic expan s ions o f functions ^'s 
for the boundary values. Many expansions of densities are found in Zolotarev ( 1986h . The method 
of the expansions may be directly applied to the projected functions ^'s. 



5 Proofs 

In this section, we give the proofs of Proposition 12.11 Theorem 12. II and Theorem 13.11 

5.1 Proof of Proposition ETTl 

For a 7^ 1, define 



1 



(2vr)'^ 70 



Then by simple notational change we have 



1 1 



1 1 



(m)"^-! e*™e-"" tan( f aes ) 

'/-^^Ad-lg™?^g-u"/cos(fQei3)e'('^/2«''s) 



du. 
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Replacing u hj u = (cos (^a^s))^^" t, we have 



2^(7ri) 



^ (cos [-aOB 



a/d 1 







vr 

cosl- 



vr 

a/d 1 
TT 



2<i(7ri)f^-i 

^ 2'^(7ri)^-i (^cos(^-o;ggJJ /i (x;a,/3ij). 

Taking the real part of this equation, we obtain the desired result when x > 0. When x < 0, replacing v 
by —V and /3 by —(3, we can obtain the result in a similar fashion. Note that ga,di~v, —(3) = ga,d{'^^ 13)- 
For a = 1, define 



1 



(2^)^^ Jo 



oo / 2 

sin(t;nH — /3?ilogn ) "dn. 



What we need is to calculate gi.di'^^P) + ilgi,d{v, P) and to consider the real part of this equation. 
Since the proof is similar to the proof of the case a 7^ 1, we omit the rest of the proof for a = 1. □ 

5.2 Proof of Theorem I01 

The function h^{x;a,P), which can be regarded as the inversion formula including the imaginary 
part, has the following representation ((2.2.20) of Zolotarev (|l986l l'): 



vr 



exp {izx + ^{z; a, — /?)) dz. 



Differentiating h^{x;a,P) n times with respect to x, we have 

1 f°° 

h^{x;a,P) = — / (iz)^ ex.p {izx + ip{z; a, —f3)) dz. 
Jo 

For the calculation of /i"'(x; a, (5) we consider the same contour as in Zolotarev ()198(tI ). which is 

z : Im {izx + Tp{z; a, —(3)) = 0, —k < arg z < — | , 



where 



-9 if ai^ 1, 
-1 if a = l, /3/-1, 
1 if a<l, /3 = -l. 



Since readers can refer to 



Lemma 2.2.3 of IZolotarev 



Zolotarevl ( 19861 ) if necessary, the details of the contour are omitted. Fr 



om 



19861 ). we only have to calculate the integration along the contour F. 



(a) (a 7^ 1, /3 7^ 1) or (a > 1, /3 = 1). Direct calculation gives 



1_ 



h^{x;a,P) = —Re / e'x.p{izx + 'il>{z;a,—P)) {iz)^dz 
vr 
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H — Im / exp {izx + 'ip{z;a, —P)) {iz)^'dz 
Jr 

-f 

+ 
1 

vr ./r 



exp {Re {izx + -i/'C-^; a, — /5))} Re 
— / exp {Re {izx + 'ilj{z;a, —(3))}lm {{iz)"dz} 
exp(-VF(¥?)) (Re + Im , 



where 



exp(-x//3)C/i(2(/?/7r;/3), a = 1. 



For details for derivation of W{ip), see p. 76 of lZolotarevI (|l986h . Replacing z by z = re"^ gives 



r e 



Since 



vr TT 
r"'(cosn((/3 + — ) + isinn(vJ + — )). 



d{r ■ sin y?) = r' sin ip dip + r cos 99 ^99, 
(i(r • cos 93) = r' cos (p dip — r sin 99 



= r' cos ip dip — r sin ip dip + i{r' sin 99 + r cos dc^. 
Combining the equations above, we have 



vr 



Re = r' sin(n + 1) + - + r cos(n + l)[ip+ — ))dip 



vr 



Im {iz)"'dz = r" ( r sin(n + 1) ( + — ^ — r' cos(n ^) ('^ ) 



From these equations and the relations of angles and the contour F on p. 76 of lZolotarevI ((1 
becomes 



Re h{x; a, j3) 



1 



TT 



exp 



_^a/(a-l)fy^(2v3/7r;^) 



vr 



xr" I r sin(n + 1) I 93 H I + r cos(n + I) [ip -\ ]] dip 



TT 



1 



exp 



-X^/'~^-^^Uaiip;9))Vn{ip)dip 



and 



Im h{x; a, /3) 



vr 



exp 



-x"/("-i)c/„(2^/vr;e) 
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xr" ( r sin(n + 1) [(p + — ^ — r' cos(n + 1) ^93 + — ) ) dip 



(b) (q = 7^ 1). Since similar results hold in this case, we omit the poof. 

(c) (a < = 1). Here we must c onsider additional contour T* = {z : Re z = 0, — r < Im z < 
as stated in p. 76 of IZolotarev (jlQSfil ). The contour T* satisfies 



Im {(izx) + ip{z; a, —1)) = 
for z gT*. The integral along T* becomes 

— / exp{izx + i/j{z;a, —l)){iz)'^dz = — / ex.p(Re {izx + 'iIj{z; a, —l))){iz)'^dz. 
vr Jp, vr 7r* 

Replacing z hy z = re~^* gives 

— j exp (Re(z2:x + ^(z; a, -1))) r"(-i)dr = -- y exp ^-ly ^-^^^ r"(ir, 

where 

(^-^^ = xrsin (^-^^ + r^cosa (^-^ + 
= —xr + r" 

for a 7^ 1 and 

= xrsin (^-^^ +rlogrsin (^-|^ + (^-^ + ^rcos (^-^^^ 
= —xr — r log r 

for a = 1. 

5.3 Proof of Proposition 13. Il 

Without loss of generality, we assume = throughout the proof. The characteristic function (j)^ 
can be written as 

<,(t) = exp(-/ |(t,s)rr(ds)+itan^ / sign(t, s) Kt, s)r r(^is) 

TTOi I \ 

-itan — y^^sign(t,s) |(t,s)| T{ds)\ 

= exp f-((j(t))" + i(a(t))" tan —/3(t)- itan— /" (t,s)r(ds)^ 
V 2 2 Jgd J 
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= exp (-(a(t))- - i{a{t)r tan ^/3(t) ((a(t))i-" - l) + i^i^t) 
= exp (-V^^^ (t)) . 

Then, the inversion formula gives 



/5(x) = (27r)- 



-i(x,t) 



exp (-V'*'^ (t)) dt 



Note that for a positive real number r > and a vector s G M'^, /ijj|(rs) = r/ijj|(s), (T(rs) = r(T(s) 
and /3{rs) = /?(s). Putting t = rs where r > and s G S'^, we obtain 



C(x) = (2vr)-'^ 



exp 



i(x, s)r - (m(s))° - i(m(s))" tan ^/3(s) ((m(s)) 



l-a 



1) 



Furthermore, replacing r by /i = r(7(s), we get 

/,^,(x) = (27r)-^/ / exp 
Js-i Jo 



^drds. 



Tra 



+ tan — /3(s) (-u 



Taking the real part of the integrand, we obtain 

■(x,s)-m1^(s) 



9a,d 



a(s) 



I3{s) (a(s))-'^ ds. 



Note that the same conclusion holds if we let a ^ 1 in the proof due to continuity at a = 1. The 
direct but tedious calculation gives the representation fifj{t) at a = 1. □ 



Acknowledgment: The authors are grateful to Prof. Makoto Maejima for very helpful comments. 
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